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ABSTRACT 


We  have  begun  to  apply  ambient  noise  surface  wave  tomography  to  broad-band  seismic  data  obtained  in  Europe, 
North  Africa,  the  Middle  East,  and  Central  Asia.  The  goal  is  to  improve  the  calibration  of  surface  wave  propagation 
in  aseismic  areas.  Ambient  noise  tomography  (ANT)  generates  surface  wave  dispersion  maps  to  periods  as  short  as 
6  across  large  areas  and,  therefore,  provides  new,  high  resolution  constraints  on  crustal  structures  homogeneously 
over  extended  regions.  Current  work  concentrates  on  the  application  of  ANT  in  Europe  where  station  coverage  is 
highest,  which  allows  improvements  in  the  method  to  be  made  rapidly  and  transported  elsewhere  later. 

In  earlier  work,  we  applied  ambient  noise  tomography  across  Europe  to  one  year  of  data  from  approximately  125 
stations  composing  the  Virtual  European  Seismic  Network  (VEBSN)  (available  from  the  ORFEUS  Data  Center)  to 
produce  Rayleigh  wave  group  velocity  maps  from  an  8-40  sec  period.  Current  work  builds  on  these  results  by  using 
longer  time  series  (to  improve  signal  quality  and  provide  measurement  uncertainties)  and  assimilating  data  from 
temporary  or  emerging  regional  networks  in  the  Mediterranean  region  to  generate  Rayleigh  wave  phase  and  group 
speed  maps  across  Europe  with  better  resolution  and  reliability.  The  study  concentrates  on  improvements  in  three 
focus  areas  where  regional  array  data  are  available:  (1)  the  western  Mediterranean,  particularly  the  Iberian 
Peninsula,  (2)  the  central  Mediterranean,  notably  Italy  and  Switzerland,  and  (3)  the  eastern  Mediterranean,  including 
Israel,  Syria,  and  Turkey.  The  western  Mediterranean  results  are  based  on  the  Spanish  National  Network  data,  which 
provides  much  tighter  constraints  on  the  crustal  structures  of  the  Iberian  Peninsula.  In  the  Central  Mediterranean,  the 
RETREAT  and  CAT-SCAN  Program  for  the  Array  Seismic  Studies  of  the  Continental  Lithosphere  (PASSCAL) 
experiments  yield  data  that  improve  resolution  in  this  area  considerably.  In  the  eastern  Mediterranean,  the  Israeli 
DESERT2000  experiment  and  the  Eastern  Turkey  PASSCAL  experiment  were  performed  simultaneously  for  an 
extended  time  period  and,  together  with  other  regional  stations,  greatly  improve  ray  coverage  in  the  region. 

We  present  Rayleigh  wave  group  and  phase  speed  measurements  with  associated  uncertainties  (where  time  series 
are  sufficiently  long)  and  dispersion  maps  in  two  of  the  three  focus  regions.  Resolution  and  the  retrieval  of 
geologically  correlated  features  on  the  dispersion  maps  are  significantly  improved  using  the  regional  data. 
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OBJECTIVE 


The  goal  of  this  research  is  to  develop  a  new  method  to  obtain  surface  wave  dispersion  measurements  based  on 
ambient  seismic  noise  and  to  produce  a  new  data  set  of  inter-station  group  and  phase  velocity  measurements.  These 
new  data  are  used  to  produce  group  and  phase  velocity  dispersion  maps,  particularly  at  periods  below  about  25  sec. 

RESEARCH  ACCOMPLISHED 


Introduction 

Surface  wave  empirical  Green  functions  (EGFs)  can  be  determined  from  cross-correlations  between  long  time 
sequences  of  ambient  seismic  noise  observed  at  different  stations.  The  feasibility  of  the  method  was  first  established 
by  experimental  (e.g.,  Weaver  and  Lobkis,  2001;  Lobkis  and  Weaver  ,  2001;  Derode  et  al.,  2003;  Larose  et  al„  2004, 
2005)  and  theoretical  (e.g.,  Snieder,  2004;  Wapenaar,  2004)  studies.  Shapiro  and  Campillo  (2004)  demonstrated  that 
Rayleigh  wave  EGFs  estimated  from  ambient  noise  possess  dispersion  characteristics  similar  to  earthquake  derived 
measurements  and  model  predictions.  Bensen  et  al.  (2007a)  describe  the  ambient  noise  data  processing  procedure  in 
detail.  The  dispersion  characteristics  of  surface  wave  EGFs  derived  from  ambient  noise  have  been  measured  and 
inverted  to  produce  surface  wave  tomography  maps  in  several  regions,  such  as  Southern  California  (Shapiro  et  al., 
2005;  Sabra  et  al.,  2005;  Gerstoft  et  al.,  2006),  the  western  US  (Moschetti  et  al.,  2007;  Lin  et  al.,  2007b),  Europe 
(Yang  et  al.,  2007),  Tibet  (Yao  et  al.,  2006),  New  Zealand  (Lin  et  al.,  2007a),  Korea  (Cho  et  al.,  2007),  Spain 
(Villasenor  et  al.,  2007)  and  elsewhere.  Ambient  noise  tomography  has  also  been  employed  on  a  continental  scale  in 
Europe  (Yang  et  al.,  2007)  and  the  US  (Bensen  et  al.,  2007b).  Most  of  these  studies  obtained  Rayleigh  wave  group 
velocity  maps  at  periods  below  a  20  sec  period,  but  longer  period  maps  at  periods  of  up  to  70  sec  (e.g.,  Yang  et  al., 
2007;  Bensen  et  al.,  2007b),  phase  velocity  maps  (e.g.,  Yao  et  al.,  2006;  Lin  et  al.,  2007b),  and  Love  wave  maps 
(e.g.,  Lin  et  al.,  2007b)  have  also  been  developed. 

Three  publications  have  resulted  from  the  research  supported  by  this  contract:  Bensen  et  al.  (2007a),  Yang  et  al. 
(2007),  and  Villasenor  et  al.,  (2007).  Bensen  et  al.  (2007a)  describes  the  data  processing  procedure  that  has  been 
developed  under  this  contract.  The  first  ambient  noise  tomography  of  the  Mediterranean  region  and  Europe  was 
performed  by  Yang  et  al.  (2007)  using  one  year  of  continuous  seismic  data  taken  predominantly  from  the  VEBSN 
obtained  from  the  Orfeus  Data  Center  in  Europe.  Attempts  to  refine  the  maps  of  Yang  et  al.  (2007)  are  based  on 
utilizing  regional  and  national  network  data  that  tend  not  to  be  available  from  the  major  international  data  centers 
such  as  the  Incorporated  Research  Institutions  for  Seismology’s  Data  Management  Center  (IRIS-DMC)  or  the 
Orfeus  Data  Center.  A  proof-of-concept  study  was  performed  using  four  months  of  data  from  the  Spanish  National 
Network  and  results  appear  in  Villasenor  et  al.  (2007). 

The  focus  of  this  paper  is  to  report  on  further  attempts  to  improve  upon  the  first  tomographic  images  provided  by 
Yang  et  al.  (2007)  by  extending  results  to  phase  velocities,  improving  resolution  by  introducing  more  seismic 
stations  predominantly  from  regional  networks  and  temporary  installations,  and  improving  the  reliability  of  the 
measurements  and  uncertainty  estimates  by  using  longer  time  series.  Because  the  available  stations  change  over 
time,  optimal  improvements  depend  on  using  data  from  several  years.  We  currently  are  working  to  process  data 
from  European  seismic  stations  acquired  from  various  sources  recorded  during  the  years  2000-2006.  Differences  in 
stations  operating  in  2000-2001  and  2005-2006  are  shown  in  Figure  1.  Regional  networks  such  as  the  Spanish 
National  Network  and  PASSCAL  experiments  are  important  resources.  The  earlier  time  periods  are  important  for 
resolution  in  the  eastern  Mediterranean  and  parts  of  the  Middle  East. 


State  of  Ambient  Noise  Data  Processing 

The  ambient  noise  data  processing  procedure  is  described  by  Bensen  et  al.  (2007a)  in  detail,  but  is  summarized  here 
briefly.  All  data  are  processed  on  a  daily  basis  and  then  are  stacked  (superposed  and  added  together)  later  to  achieve 
cross-correlations  over  arbitrarily  long  time  series  currently  ranging  up  to  two  years.  The  mean,  trend,  and 
instalment  response  of  the  daily  three-component  (E,  N,  Z)  seismograms  are  removed  and  the  resulting  waveforms 
are  band-pass  filtered  between  periods  of  5  sec  and  100  sec.  We  now  compute  the  full  cross-correlation  tensor  and 
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not  just  vertical-vertical  cross-correlations.  To  speed  up  the  process,  we  do  not  rotate  the  components  into  the  radial 
(R)  and  transverse  (T)  directions  for  each  station-pair  until  the  component  cross-correlations  (E-E,  E-N,  N-N,  N-E) 
are  performed.  Earthquake  signals  and  instrumental  irregularities  are  then  removed  by  temporal  normalization.  In 
order  to  postpone  the  component  rotation  until  after  cross-correlation,  the  East  and  North  components  are  temporally 
normalized  together.  To  achieve  this,  both  components  are  band-pass  fdtered  between  15  sec  and  50  sec,  a  band  that 
contains  the  most  energetic  surface  wave  signals  from  earthquakes.  For  each  time  point,  the  mean  of  the  absolute 
value  of  each  seismogram  is  computed  in  the  128  sec  window  centered  on  that  point.  The  values  of  the  East  and 
North  components  are  compared,  and  the  larger  is  used  to  define  the  inverse  weight  for  that  time  point.  That  weight 
is  then  applied  to  both  the  North  and  East  component  time-series,  band-passed  between  5  sec  to  100  sec.  This 
process  effectively  suppresses  earthquake  signals  and  maintains  the  linearity  of  the  rotation  operator.  After  temporal 
normalization,  the  signals  are  whitened  in  frequency.  Before  whitening,  ambient  noise  is  most  energetic  in  the 
microseism  bands  below  20  sec  period.  Frequency  whitening  is  carried  out  to  broaden  the  period  band  of  the 
dispersion  measurement.  Again,  to  maintain  the  linearity  of  the  rotation  operator,  the  East  and  North  signals  are 
whitened  together.  Because  the  spectra  of  both  components  are  similar,  on  average,  we  do  this  by  simply  weighting 
the  East  and  North  signals  in  the  frequency  domain  by  the  inverse  of  the  smoothed  East  spectrum.  Other  methods, 
such  as  weighting  by  the  mean  of  the  two  spectra  or  their  product,  produce  similar  results.  We  then  stack  all 
available  daily  cross-correlations  for  each  component  and  station-pair  in  order  to  enhance  the  signal-to-noise  ratio 
(SNR).  Because  all  of  these  processes  are  linear  in  the  rotation  operator,  the  transverse-transverse,  transverse-radial, 
radial-radial,  and  radial-transverse  cross-correlations  between  each  station-pair  can  be  calculated  by  a  linear 
combination  of  those  four  components.  The  resulting  cross-correlations  can  be  converted  to  Empirical  Green’s 
functions  by  introducing  an  additive  phase  factor  (Lin  et  al.,  2007b). 

Examples  of  vertical  component  cross-correlation  record  sections  are  shown  in  Figure  2.  Figure  2a  is  based  on  four 
months  of  ambient  noise  data  from  the  Spanish  National  Network.  Figure  2b  is  from  12-month  time  series  centered 
on  a  station  in  central  Italy.  The  arrivals  near  zero  cross-correlation  lag  are  not  interpreted  in  this  study,  and  their 
origin  remains  unknown.  We  believe  that  they  result  from  the  conversion  of  ocean  waves  to  mantle  propagating 
P-waves  in  the  European  Arctic.  They  are  observed  elsewhere  across  Eurasia,  and  understanding  their  generation 
and  nature  requires  further  effort.  They  do  not,  however,  interfere  with  the  fidelity  of  the  surface  wave  dispersion 
measurements. 

Yang  et  al.  (2007)  processed  one  year  of  data  from  2004  taken  from  about  120  stations  available  from  the  Orfeus 
Data  Center  in  Europe.  Station  coverage  was  similar  to  the  red  symbols  shown  in  Figure  la. 

Improvements  over  the  results  of  Yang  et  al.  (2007)  derive  from  three  factors.  First,  we  will  use  longer  time  series, 
two  years  of  data,  which  will  allow  the  computation  of  better  measurement  uncertainties.  Uncertainties  are 
computed  from  repeating  the  measurements  in  different  seasons.  Second,  there  are  data  resources  in  Europe,  North 
Africa,  and  the  Middle  East  that  were  not  used  by  Yang  et  al.  These  include  national  networks  such  as  the  Spanish 
National  Network  and  the  Israeli  DESERT2000  experiment,  PASSCAL  experiments  such  as  the  RETREAT  and 
CAT-SCAN  experiments  in  Italy  and  Eastern  Turkey  and  Northern  Anatolian  Fault  (ANF)  experiments  in  Turkey, 
as  well  as  other  temporary  deployments  such  as  the  Mid-Sea  Experiment  and  the  Carpathian  Basin  Project  (CBP)  in 
Austria  and  Hungary.  We  are  now  in  the  process  of  acquiring  and  beginning  to  process  these  data.  Focus  is  being 
placed  on  three  regions:  (1)  the  western  Mediterranean,  (2)  the  central  Mediterranean,  and  (3)  the  eastern 
Mediterranean.  Third,  we  have  extended  the  measurements  to  include  phase  velocities.  Lin  et  al.  (2007b)  show  that 
phase  velocity  measurements  are  more  repeatable  than  group  velocities  and,  thus,  have  lower  uncertainties. 

Western  Mediterranean:  Villasenor  et  al.  (2007)  computed  cross-correlations  between  stations  within  the  Spanish 
National  Network  (SNN)  during  July,  August,  September,  and  October  of  2005.  Initially,  only  group  velocity 
tomography  maps  were  obtained,  but  we  have  subsequently  measured  phase  velocities  from  these  data.  Additionally, 
group  and  phase  speed  measurements  were  obtained  from  cross-correlations  between  VEBSN  and  GEOFON 
stations  on  the  Iberian  Peninsula  (available  from  Orfeus  and  IRIS,  respectively)  operating  from  April  2004  to  March 
2005. 

Central  Mediterranean:  We  have  augmented  the  data  processing  performed  by  Yang  et  al.  (2007)  by  computing 
cross-correlations  and  obtaining  group  and  phase  velocities  from  one  year  of  continuous  data  (April  2004-March 
2005)  using  the  RETREAT  and  CAT/SCAN  PASSCAL  experiments  in  Northern  and  Southern  Italy,  the  VEBSN, 
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and  MEDNET  station  data  available  from  IRIS.  The  addition  of  data  from  the  two  PASSCAL  experiments  in  Italy 
significantly  improves  upon  the  initial  results  found  by  Yang  et  al.  (2007). 

The  Eastern  Mediterranean/Middle  East:  We  are  currently  processing  data  from  stations  that  operated  during 
2000  and  2001.  This  data  set  includes  stations  from  the  Eastern  Turkey  (PASSCAL)  Seismic  Experiment,  the 
DESERT2000  experiment  in  Israel,  and  stations  available  from  IRIS  in  Western  Turkey  and  Crete. 


Current  Status:  Example  Group  and  Phase  Velocity  Maps 

Figures  3  and  4  present  examples  of  the  current  state  of  the  group  and  phase  velocity  maps  across  the  region  of 
study.  The  tomographic  method  of  Barmin  et  al.  (2001)  has  been  used.  One  aspect  of  this  method  is  that  outside  the 
region  of  high  path  coverage,  the  maps  revert  to  the  average  of  the  measurements.  Thus,  the  periphery  of  all  of  the 
maps  are  white  because  in  these  regions  there  are  few  ray  paths.  Tomography  is  performed  on  a  1  deg  x  1  deg  grid 
over  the  region  of  study.  Figure  5  is  a  map  of  sedimentary  and  crustal  thicknesses  taken  from  CRUST2.0,  for 
comparison. 

In  Figure  3,  low  velocities  on  the  10  sec  maps  reveal  the  location  of  major  sedimentary  basins.  The  Po  Basin  in 
Italy,  the  North  Sea  Basin,  the  Silesian  Basin  in  northern  Germany  and  Poland,  the  Tyrrhenian  and  Adriatic  Basins 
surrounding  Italy,  and  the  Valencia  Trough  east  of  Spain  are  some  examples.  The  Valencia  Trough  and  the  Sub- 
Pyrenees  Basin  are  seen  better  in  Figure  4.  At  a  25  sec  period,  the  maps  reveal  information  about  vertically 
integrated  crustal  shear  velocities  and  crustal  thickness  (the  latter  particularly  in  phase  velocities  that  sample  deeper 
than  group  velocities).  Low  velocities  are  associated  with  thicker  crust.  Thus,  there  are  low  velocities  beneath  the 
Betic  Mountains  in  southern  Spain,  the  Alps,  the  Carpathians,  the  Dinarides  in  the  Balkans,  and  the  Apennines  along 
the  spine  of  Italy.  The  Hellenic  Arc  also  appears  as  low  velocities,  but  probably  not  because  of  particularly  thick 
crust  there. 


Expected  Continued  Improvements 

Figure  6  shows  how  path  coverage  has  improved  with  the  addition  of  measurements  obtained  in  Spain  using  the 
Spanish  National  Seismic  Network  and  data  from  two  PASSCAL  experiments  in  Italy  (RETREAT,  CAT/SCAN). 
Nevertheless,  the  high-resolution  (non-white)  regions  of  the  maps  in  Figures  3  and  4  are  mainly  within  continental 
Europe.  Our  goal  is  to  improve  data  coverage  in  the  Mediterranean  Sea,  North  Africa,  and  parts  of  the  Middle 
East.  We  plan  to  continue  to  improve  the  data  set  in  each  of  the  three  regions  by  increasing  the  number  of  stations, 
and  increasing  the  time  duration  of  ambient  noise  analysis. 

(1)  We  will  improve  upon  the  initial  results  of  Villasenor  et  al.  (2007)  by  increasing  the  time  duration  of  analysis  of 
the  Spanish  National  Network  from  four  months  to  two  years  (2005  and  2006).  Additionally,  we  will  include  in  our 
analysis  data  from  VEBSN,  IRIS,  MEDNET  and  other  stations  operating  in  Europe  during  this  time  period.  By  tying 
the  Spanish  National  Network  data  to  VEBSN  and  other  stations,  resolution  in  the  western  Mediterranean  will  be 
improved  dramatically. 

(2)  To  improve  ray  path  coverage  and  resolution  in  the  central  Mediterranean  we  will  compute  cross-correlations 
between  stations  operating  during  2005  and  2006,  which  includes  the  RETREAT  and  CAT/SCAN  PASSCAL 
experiments  in  Italy  and  the  Spanish  National  Network,  as  well  as  VEBSN  and  other  stations  located  in  this  region 
(Fig  lb).  These  data  will  result  in  additional  ray  coverage  across  the  Mediterranean  Sea  between  Italy  and  the 
Iberian  Peninsula  and  between  Italy  and  the  Hellenic  Arc. 

(3)  Cross-correlations  between  stations  operating  during  the  years  2000  and  2001  from  the  Eastern  Turkey  Seismic 
Experiment,  the  DESERT2000  experiment  in  Israel,  and  stations  available  from  IRIS  in  Western  Turkey  and  Crete 
will  yield  improved  resolution  across  Turkey,  Syria,  Israel,  and  the  Eastern  Mediterranean. 
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CONCLUSIONS 


Three  studies  supported  by  this  contract  have  already  been  published.  Bensen  et  al.  (2007a)  described  ambient  noise 
data  processing  in  detail  and  Yang  et  al.  (2007)  performed  ambient  noise  tomography  for  Rayleigh  wave  group 
velocities  using  observatory  data  in  Europe  that  was  performed  between  10  sec  and  40  sec.  In  addition,  assimilation 
of  regional  data  by  Villasenor  et  al.  (2007),  notably  using  the  Spanish  National  Network,  established  the  ability  of 
such  data  to  improve  resolution  at  short  periods  (6-20  sec)  when  used  in  concert  with  observatory  data.  Further  data 
processing  using  PASSCAL  experimental  data  in  Italy  has  established  similar  results  in  the  central  Mediterranean 
region,  and  the  data  processing  methodology  has  been  extended  to  include  phase  velocities.  In  the  near  future,  more 
regional  data  will  be  assimilated  into  the  inversion  to  improve  data  coverage  and  resolution  across  much  of  the 
Mediterrean  Sea,  in  particular  in  the  eastern  Mediterrean,  including  parts  of  the  Middle  East  and  North  Africa. 
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2000-2001 


2005  -  2006 


Figure  1.  Contrast  in  station  coverage  in  the  years  2000-2001  versus  2005-2006.  The  station  color  code  is  as 
follows:  BLUE — regional  networks  (e.g.,  Spanish  National  Network,  the  Israeli  DESERT2000 
experiment);  RED — other  data  available  through  IRIS;  YELLOW — VEBSN  data  available  from 
the  Orfeus  Data  Center  in  Europe;  PURPLE — PASSCAL  experiments  or  other  European 
PASSCAL-like  deployments  (such  as  the  MIDSEA  experiment);  ORANGE — SeisUK  stations. 
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Figure  2.  Record  sections,  (a)  From  the  Spanish  National  Network.  Cross-correlations  are  from  4  months  of 
data,  (b)  From  stations  around  Italy.  Cross-correlations  are  from  one  year  of  data. 
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Figure  3.  Examples  of  the  current  status  of  Rayleigh  wave  group  and  phase  velocity  maps  across  Europe. 

Phase  velocities  appear  in  the  left  column  and  group  velocities  in  the  right  column.  Maps  at  a  10  sec 
period  are  in  the  top  row  and  at  a  25  sec  period  in  the  bottom  row. 
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Figure  4.  Rayeigh  wave  10  sec  group  velocity  map  across  Spain  presented  with  a  tectonic  map  for 

comparison.  Low  velocities  are  associated  with  the  sub-Pyrenees  Basin  and  the  Valencia  Trough 
off-shore.  High  velocities  in  the  west  are  associated  with  the  Paleozoic  Iberian  Massif  and  the 
Iberian  Range  is  also  correlated  with  high  velocities  at  this  period. 
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Figure  5.  Sediment  thickness  and  crustal  thickness  presented  for  comparison  with  Figure  4.  Values  taken 
from  CRUST2.0. 
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Figure  6.  Path  density,  (a)  From  the  study  of  Yang  et  al.  (2007).  (b)  Current  state  of  path  density,  showing 
improvements  in  Spain,  central  Europe,  Italy  and  the  central  Mediterranean  from  recent  data 
processing.  Path  density  is  defined  as  the  number  of  paths  intersecting  a  local  2  deg  X  2  deg  cell. 
Addition  of  more  data  from  other  regional  networks  will  improve  path  density  across  the  entire 
Mediterranean,  particularly  in  the  eastern  Mediterranean  and  parts  of  the  Middle  East. 
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